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Abstract 

In this work we address the open problem of high Reynolds number limit in hydrodynamic 
turbulence, which we modify by considering a vanishing random (instead of deterministic) viscosity. 
In this formulation, a small-scale noise propagates to large scales in an inverse cascade, which can 
be described using qualitative arguments of the Kolmogorov-Obukhov theory. We conjecture that 
the limit of the resulting probability distribution exists as Re —>■ oo, and the limiting flow at 
finite time remains stochastic even if forcing, initial and boundary conditions are deterministic. 
This conjecture is confirmed numerically for the Sabra model of turbulence, where the solution is 
deterministic before and random immediately after a blowup. Then, we derive a purely inviscid 
problem formulation with a stochastic boundary condition imposed in the inertial interval. 
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I. INTRODUCTION 


The phenomenon of turbulence remains a major challenge in incompressible fluid dynam¬ 
ics. In the idealized formulation, when forcing is applied at a large scale L, the developed 
turbulence is considered to be statistically stationary, isotropic, homogeneous, and universal 
at small scales. The contemporary understanding of this phenomenon largely relies on the 
Kolmogorov-Obukhov theory Bca. which conjectures that the mean energy flux £ from large 
to small scales together with the viscosity v define all statistical properties of the small-scale 
dynamics. The problem of developed turbulence contains a natural dimensionless small 
parameter, the inverse Reynolds number. Since the Reynolds number achieves very large 
values in applications, e.g.. Re ~ 10® in atmospheric flows, the theory accessing the limiting 
description as Re —)■ oo would be the most appropriate. 

The problem of describing the limit of large Reynolds numbers remains open. Mathe¬ 
matically, weak solutions of the Euler equations should be associated to the inviscid limit 
of the incompressible 3D Navier-Stokes equations; however, extra admissibility conditions 
are required due to unphysical solutions |3H5]. A different approach involves a probabilistic 
description. The necessity of a stochastic approach for turbulent flows is widely accepted j6] 
and justihed, e.g., by the exponential path separation controlled by positive Lyapunov ex¬ 
ponents. However, the dynamical system viewpoint fails in the limit Re —>■ oo, since the 
Lyapunov exponents depend strongly on Re. The stochastic formulation is traditionally 
achieved, e.g., by considering random large-scale forcing. Though the relevance of random 
forcing as a physical model can be questioned, it is a convenient theoretical assumption. On 
the other hand, the recent ideas of spontaneous stochasticity developed for Lagrangian tra¬ 
jectories in fully developed turbulent flows [THIQ] suggest the randomness as a consequence 
of non-uniqueness of solutions. Hence, the stochasticity may arise naturally in the equa¬ 
tions, which look fully deterministic (deterministic forcing, initial and boundary conditions) 
as soon as the limit Re —)■ cx) is considered. 

This work is aimed to extend the ideas of spontaneous stochasticity, at the same time 
following the earlier studies on the influence of a small-scale noise [nHH]. We provide a 
numerical evidence and a qualitative theoretical reasoning suggesting that the turbulent flow 
in the limit Re —?• oo is an intrinsically stochastic process well defined within a deterministic 
formulation, i.e., for deterministic forcing, initial and boundary condition. Such a descrip- 
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tion is achieved by introducing an infinitesimal (vanishing as Re —>■ oo) random perturbation 
of the viscosity, acting as a regularizing term in the governing equations. We show that this 
perturbation acts at the Kolmogorov viscous scale and induces an inverse cascade resulting 
in stochastic mixing at larger and larger scales m- This process can be described qualita¬ 
tively using arguments of the Kolmogorov-Obukhov theory. As a result, the stochastic flow 
components persist at all scales in the limit of vanishing random viscosity. 

The conjecture that we call the stochastic anomaly and verify numerically for the Sabra 
model of turbulence suggests the existence of a probability distribution in the limit Re —?• oo, 
which is independent of the small-scale regularizing noise. Our numerical limiting solution 
is deterministic at times before the blowup, and it becomes random spontaneously at later 
times. On one hand, this means that the flow cannot be predicted exactly at a finite time, 
as opposed to the chaotic motion where a finite-time uncertainty can always be removed 
for sufficiently accurate initial conditions. On the order hand and more importantly, this 
shows that the limiting (Re —>■ cx)) probability distribution is well-defined and, thus, can be 
predicted. We demonstrate theoretically and numerically that the limiting probability dis¬ 
tribution satisfies the equations of fully inviscid dynamics with a proper stochastic boundary 
condition imposed in the inertial interval. 

In Section [TT[ we give a qualitative description of the inverse cascade of stochasticity, 
based on Kolmogorov arguments. Section uni formulates the hypothesis of the stochastic 


anomaly in the limit of large Reynolds numbers, and Section IV verifies this hypothesis 
numerically for the Sabra model of turbulence. Section |V] describes main features of the 
suggested stochastic description, as compared to the phenomena of deterministic chaos and 


spontaneous stochasticity of Lagrangian trajectories. Section VI describes the small-scale 
statistics, which is used in Section VII| for implementing the purely inviscid formulation with 
a stochastic closure at small scales. Section VIII| provides numerical tests of this probabilistic 
formulation. We finish with some conclusions. 


II. THE INVERSE CASCADE OF STOCHASTICITY 

The importance of small-scale noise for describing a turbulent flow was recognized for 
long time, see e.g. [TT1 - IT3] . Due to large separation between the scale, where the noise is 
“injected”, and its eventual observation at large scales of the flow, this becomes a cascade 
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phenomenon. As a result, the noise amplihcation is governed by a power-law, as opposed to 
the exponential separation of trajectories in hnite-dimensional deterministic chaos. In this 
section, we provide a qualitative argument based on the Kolmogorov approach describing 
the propagation of stochastic perturbations from small to large scales. 

The Kolmogorov-Obukhov theory BE] conjectures that the mean energy flux £ from 
large to small scales together with the viscosity v dehne all statistical properties of the 
small-scale dynamics. At scales of the inertial interval, t] i L, where 

^ = Z.3/4^-1/4 

is the Kolmogorov length, both forcing and viscosity can be neglected. The dimensional 
analysis yields power-laws for the moments of velocity fluctuations, (|(5'y|^) ~ see 

e.g. [IS]- These scaling laws are only approximate due to the existence of anomalous cor¬ 
rections in the exponents (Bj, which still lack the comprehensive theoretical explanation. 
At smaller scales, i rj, the flow is dominated by viscosity. Choosing the large-scale 
speed as K = the Reynolds number is related to the inertial interval span as 

Re = VL/u = 

Assuming locality of nonlinear interactions, the energy produced by forcing at scale L 
is transferred successively to smaller and smaller scales, until it dissipates at the smallest 
viscous scales. In the inertial range, a turnover time at scale I is given by r i/ 3£2/3^ and 
it is a characteristic time for the energy transfer to the smaller scale £/2. Hence, a total time 
for the energy transport from the forcing to viscous scales is given by a convergent geometric 
series, independent of viscosity. This observation constitutes the dissipation anomaly due to 
Onsager HSlEj: the energy flux to small scales remains hnite for vanishing viscosity, i.e., 
in the limit of high Reynolds numbers. Re —)■ oo. 

A similar phenomenon occurs for a noise, but with an opposite direction from small to 
large scales nulla]. There are various physical reasons for introducing a noise in turbu¬ 
lent flow at small scales, e.g., thermodynamic molecular motion or parameter uncertainties. 
Apparently, an effect of this noise at large scales is independent on its particular form and 
origin. In this work, we model the noise by considering an uncertainty for the viscosity 
parameter (1 -|- x)v, where x > —1 is a random number that describes a relative viscosity 
deviation. Such an uncertainty may result from various physical properties: humidity fluctu¬ 
ations in air or temperature variations in water, etc. The measurement and numerical errors 
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also serve as a justification. This viscosity perturbation yields an additional viscous term 
xz/Av in the governing Navier-Stokes equations. In the inertial interval, one can use the 
Kolmogorov estimate Sv ~ for the speed variation and for the Laplace operator, 

which leads to the relation 


xz/Av ~ 


( 2 ) 


where we used Eq. (|^ and denoted by the standard deviation of the relative viscosity 
perturbation x. Similarly, the quadratic convective term (v ■ V)v has the scaling 
This shows that, for cxa; ^ 1, the random viscous term can be neglected at scales of the inertial 
interval, £ S> ? 7 , just like one neglects the deterministic viscosity. Thus, the random viscous 
term becomes important only at the Kolmogorov scale, i ^ rj. This can be interpreted as 
the “injection” of a stochastic component to the flow at a small-scale end of the inertial 
interval. 

The effect of random viscosity perturbation on the dynamics at Kolmogorov scale can be 
estimated from the Newton’s law, i.e., comparing a random viscous force in Eq. ^ at i = rj 
with the acceleration term D\/Dt. The stochastic velocity component grows in magnitude 
to the mean value (|5n|) ~ the characteristic time 


tx 


< 

rsj 


(|hn|) 


i2/3 


a. 


,77-1/352/3 


( 3 ) 


Considering the large-scale characteristic time as T = LjV = e i/^L^/^ yields 

1 


T ~ a. \l) 


a. 




( 4 ) 


This means that the time of development of stochastic dynamics at the Kolmogorov scale is 
negligible for the large-scale flow if 


CTt, 3> 


V^Re 


( 5 ) 


In the limit Re —?• 00 , Eq. ([^ ensures that even a vanishingly small relative perturbation 
of the viscosity is sufficient for our further arguments. We note that Eq. (|^ yields a rather 
mild estimate based on the linear growth of disturbances. Similar arguments taking into 
account the exponential path separation yield much smaller values of tx na. 

Injection of a random flow component at small scales of the inertial interval leads to 
its turbulent transport towards the largest scale L due to nonlinear interaction mm, 
i.e., to an inverse cascade of stochasticity. In this process, the time required for developing 
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the stochastic velocity component at scale 21 due to analogous component at scale I is 
comparable to the turnover time r ~ which rapidly increases with the transition 

to larger and larger scales. Because of the power-law relation for r, the total time of the 
transition through the inertial interval is given by a convergent geometric series, i.e., it 
remains hnite in the limit Re —)■ cxd and does not dependent on viscosity na. Also, due to 
the increase of turnover times for larger scales, we expect that the fast small-scale stochastic 
dynamics is self-averaged on top of the slow large-scale flow, losing the dependence on a 
probability distribution of the random variable x. 


III. STOCHASTIC ANOMALY 

Reynolds numbers appearing in applications of the developed turbulence attain very 
large values. For example, in typical atmospheric flows. Re ~ 10® results in the huge scale 
separation between the viscous and large-scale motion, L/rj Re®/^ 10®. Since the large- 

scale flow is usually of primary interest, it is reasonable to consider the limit of high Reynolds 
numbers. In this section we formulate a hypothesis of the stochastic anomaly, which is a 
consequence of the inverse cascade of stochasticity considered in the limit Re —?• oo. This 
hypothesis will be conhrmed numerically for the Sabra shell model of turbulence in the 
following sections. 

The two properties of the inverse cascade of stochasticity are important for taking the 
limit Re —?• oo. First, we argued that specihc properties of the small-scale noise are “for¬ 
gotten” at large scales, due to an increasing turnover time. Second, we saw that the arrival 
time of a stochastic perturbation to a specihc scale has a hnite limit for high Re. This 
suggests that the limit Re —)• oo should be dehned in terms of probability distributions for 
the velocity held. In this dehnition, the random viscosity plays a role of a regularization 
term, similarly to the role of a deterministic viscosity in the Burgers equation leading to a 
shock wave in the inviscid limit. 

Therefore, our conjecture states that, in the limit of high Reynolds numbers formulated 
for a vanishing random viscosity, a limiting how is described by a time-dependent probabil¬ 
ity distribution determined for given deterministic forcing, initial and boundary conditions. 
Since this limiting solution is stochastic, despite all the conditions of the how are determin¬ 
istic, we call it the stochastic anomaly. 
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IV. NUMERICAL SIMULATIONS FOR SABRA MODEL 


Direct numerical simulations of the incompressible 3D Navier-Stokes equations, which 
accurately resolve all scales involved in the dynamics, are limited to rather low Reynolds 
numbers due to computer limitations. Still, condition ([^ shows that verihcation of the 
stochastic anomaly is possible at the boundary of its validity. For example. Re = 100 
requires ^ 0.1, which is satished for ^ 1 when the random viscosity perturbation xu 
remains small. However, computing the probability distribution requires a large number of 
simulations, each for a different random value of x, leading to extremelly high requirements 
for computation resources. 

Instead, in this work we perform the detailed verihcation using the Sabra model of turbu¬ 
lence [TS], which is a natural playground for the studies on developed turbulence. This model 
allows accurate numerical simulations for very high Reynolds numbers and possesses non¬ 
trivial properties of the Kolmogorov-Obukhov theory: energy cascade, dissipation anomaly, 
scaling power laws for velocity moments and anomalous corrections close to the ones for the 
Navier-Stokes equations. The Sabra model is obtained by reducing dynamics to a discrete 
sequence of shells |k| = in the Fourier space for the geometric progression of wavenumbers 
kn = fcoA"", n = 1,2,3,... (we use = 1 and A = 2). The turbulent “how” is described 
by complex velocity variables Un{t), which mimic the velocities at the corresponding shells. 
Thus, the velocity characterizes the scale i ~ k~^. The model equations are 


'kn 


ikn ^^'^n+2'^n+l 2^^n+l'^n—l ^'^n—l'^n—2 

-(1 + +/„, n = l,2,3, 


( 6 ) 


where /„ is the forcing and [1 + x)v is the viscosity with a random perturbation; the 
boundary conditions are chosen as uq = U-i = 0. Eq. 0 has the inviscid invariants: energy 
E = \un\‘^ and helicity H = ^(— 

In our simulations, we choose the large-scale deterministic non-stationary forcing /2 = 
and /a = 1 — f and initial conditions Ui = 1/2 + i, U 2 = 1 dX t = Q (all other components 
are zero) in the model with 40 shells. In this case, the characteristic values L, V and T 
are of order 1 (a.u.), and the Reynolds number is defined as Re = 1/u. The values up to 
Re = 10^° are considered. In each test, numerical simulations are performed for 10^ different 
values of the viscosity (1 -|- x)^, where x (hxed for each simulation) is taken as a random 
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variable uniformly distributed in the interval |x| < a/ScTj. with ax = 0.3. Additional tests 
with different initial conditions, forcing and distributions for x led to the same conclusions. 

Fig. (the hrst three columns) presents the numerical probability densities for the large- 
scale velocity U 2 {t) at different times and for different Reynolds numbers. One can clearly 
see the stochastic form of large-scale dynamic, which converges with the increase of Re. 
This conhrms our hypothesis of the stochastic anomaly suggesting that, for high Reynolds 
numbers, the flow remains stochastic with a well-dehned limiting probability distribution. 
Note, however, that the convergence is rather slow: one can clearly recognize the limiting 
pattern for Re = 10^, while the distortion is still rather strong for Re = 10®. In Fig. |^we 
compare mean values and standard deviations of the large-scale speeds for Re = 10^ and 
10^°. We stress again that the presented probability distributions correspond to specihc 
deterministic initial and forcing conditions. 

Fig. I^b) shows that the stochastic components of the flow do not exist up to the time 
th ~ 0.42. This time represents the blowup in the inviscid Sabra model, when the enstrophy 
grows to inhnity and velocities at all scales get a power-law excitation (THl EU] . 
Such behavior fully agrees with the regularity result proving uniqueness of the solution 
for the inviscid Sabra model at times before the blowup. Nonzero standard deviations for 
t > tf, identify the blowup as a source of the spontaneous stochasticity: the solution is 
deterministic for t < tf, and random for t > th¬ 
in this section we presented the large-scale shell speeds, whose behavior is the most 
important from practical point of view. However, the nature of turbulent flow is in small- 
scale fluctuations [SI US]- We return to the study of the probability distribution for small- 
scale velocities in Section lYD 

V. STOCHASTIC ANOMALY VS. DETERMINISTIC CHAOS AND SPONTA¬ 
NEOUS STOCHASTICITY OF PARTICLES 

It is common to think that the developed turbulence cannot be described in terms of the 
theory of deterministic chaos. The spontaneous stochasticity of Lagrangian trajectories [Zh 
[To] is one of the properties that highlights the conceptual difference. In this section we 
discuss these issues from the point of view of the stochastic anomaly. 

In hnite-dimensional deterministic chaos, the concept of large-time unpredictability is 
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FIG. 1. Probability density, from blue (zero probability) to red (maximum probability) of the 
velocity U 2 {t) on complex plane at different times (rows) and Reynolds numbers (columns), obtained 
for deterministic initial conditions and forcing. The last column is obtained using the inviscid 


stochastic formulation of Section VII with only 7 shells. 


based on exponential separation of solutions, which start at close initial points. A rate of 
separation is determined by positive Lyapunov exponents. However, the system dynamics is 
deterministic at every finite time: sufficiently (exponentially) close initial conditions yield a 
hnite-time prediction with an arbitrary accuracy. The probabilistic description is introduced 
by considering a chaotic attractor in the limit t —)■ oo, i.e., an invariant probability measure. 
This stochastic approach is applicable at sufficiently large times describing a statistical 
equilibrium. Such properties clearly distinguish the deterministic chaos from the stochastic 
anomaly: The limiting (Re —)■ cxd) turbulent flow is truly stochastic at hnite times, even 
if the initial conditions are fully deterministic. Furthermore, its probability distribution is 
time-dependent, not stationary. 


9 















(a) (b) 




FIG. 2. (a) Mean values and (b) standard deviations for real parts of the large-scale shell speeds 
ui, U 2 and U 3 as functions of time. The speeds become stochastic after the blowup at t = 0.42. 
Solid black curves correspond to Re = and dashed red curves to Re = 10^. 

In order to see how this is possible, let us consider a simple ordinary differential equation 

(7) 


'f = 


which mimics a particle position r at time t with the initial condition r(0) = 0. Here the 
velocity v = is chosen such that it satishes the Kolmogorov scaling. There is a family 
of solutions 

0 t <ts] 


r{t) = 


3 


3/2 


( 8 ) 


, t > tg, 


where tg is an arbitrary parameter denoting a spontaneous time, when the particle starts 
moving, see Fig. This example shows the non-uniqueness of the trajectories, inherent in 
Kolmogorov scaling laws. For Eq. Q, this non-uniqueness is the well-known fact in dif¬ 
ferential equations, because the right-hand side if not Lipschitz continuous. For turbulent 
flows, these ideas appeared and were further elaborated for fluid particle trajectories assum¬ 
ing given (hxed or stochastic) rough velocity helds, see e.g. [TUTU]. As one can see from 
Eq. (|^, a separation between two solutions with close initial conditions grows as power-law, 
not exponentially. Moreover, if one of the initial conditions is at the origin, different sepa¬ 
rations can be achieved at a given time t by choosing different parameters tg. This shows 
how the stochasticity emerges instantaneously due to the non-uniqueness, as opposed to 
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FIG. 3. Non-unique solutions of Eq. Q with the initial condition r(0) = 0 (solid lines) span a 
gray region, which grows as a power-law with time. The dashed line shows a solution with a small 
non-zero initial condition. 


the exponential path separation in the deterministic chaos. Note also that the spontaneous 
stochasticity does not require the system to be chaotic, as follows from the above example 
and can be observed in the Gledzer shell model of turbulence [22] . 


Both in Eq. ([^ and in the theory describing non-unique particle trajectories [THin], a 
singularity is introduced explicitly in the governing equations. A more sophisticated process 
drives the stochastic anomaly as described in Section IV Here, equations ([^ do not feature 
any singularity in the right-hand side. Instead, a singularity leading to the non-uniqueness 
appears in the solution itself at the hnite-time blowup. 


Though the unpredictability of spontaneously stochastic turbulent flow is qualitatively 
different from the unpredictability in deterministic chaos, both lead to similar practical 
conclusions on essential limitations for hnite-time predictions. The stochastic anomaly, 
however, comes along with a solution to this problem: it suggests that the regular probability 
distribution exists in the limit Re —)■ oo. In this sense, the spontaneous stochasticity is 
a property already inherent in the inviscid how equations, i.e., inviscid Sabra model or 
incompressible Euler equations. This makes the limiting probability distribution a true 
physical solution of a “deterministic” inviscid problem that can be computed as a function 
of time and, thus, accurately predicted. In the following sections we suggest how this can 
be done in the framework of the Sabra shell model. 
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FIG. 4. (a) Standard deviations of real parts of shell velocities (logarithmic vertical scale); at last 
time, the shell numbers n = 1, 2,... are ordered from top to bottom, (b) Instantaneous average of 


the dissipation rate D = 2(1 + x)h'Y^ /c^|unp. (c) Velocity moments \un{t)\P at times t = I, 1.5, 2 
(solid, dashed and dash-dot lines) for p = 1, 2,..., 5 (black, red, blue, etc.). The graphs are based 


on the solutions with Re = see Section IV 


VI. SMALL-SCALE BEHAVIOR 


In this section, we focus on the time-dependent statistics at small-scales, i.e., in the 
inertial interval, for the simulations of Section IV Fig. |^a) shows the standard deviations 
for real parts of the shell speeds. One can see that all these quantities become non-zero 
after the blowup and change continuously with time (a small noise on top of the curves has 
a numerical origin due to finite sampling). The forcing range of shells, n ^ 4 at later times, 
can be clearly distinguished. In the inertial interval corresponding roughly to the shells 
n ^ 5, the standard deviations accompany a decrease of the mean dissipation rate in time 
followed by a small increase, as one can see in Fig. |^b). Recall that our simulations use 
non-st at ionary deterministic forcing and, hence, no stationarity is expected in Fig. |^a,b). 

For stationary turbulence, the scaling of time-averaged velocity moments {\un\^) ~ kn'^^ 
is known to be anomalous, deviating from the Kolmogorov prediction (p = p/3 [23]. This 
reflects the intermittency phenomenon, characterized by strong short-time bursts of shell 
speeds separated by regions of slow behavior. In our simulations, the probability distribution 
is non-st at ionary, changing with time at all scales. Considering the averaging at a hxed 
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time denoted with a bar (instead of the common time averaging), we obtain an important, 
though expected, observation that the scaling laws are valid instantaneously with the same 
exponents, see Fig. |^c). In the figure, the graphs at times t = 1, 1.5 and 2 are presented 
in the logarithmic scale and have the same slope for each p (vertical shifts are used for 
better comparison). These slopes shown by gray lines are given by the anomalous exponents 
Cl = 0.39, C 2 = 0.72, Cs = 1, C 4 = 1-26, Cs = 1-48 of the stationary turbulence [18] . 

Our numerical observations suggest that the stationarity assumption for the description 
of the inertial interval can be relaxed, if one averages over the instantaneous (spontaneously 
stochastic) distribution, instead of performing the time averaging. Furthermore, this dis¬ 
tribution depends slowly with time (except in a vicinity of the blowup). Thus, the in- 
termittency remains a property of a single realization of the stochastic process, while the 
time-dependence of the probability distribution at all scales is characterized by the large- 
scale time T determined by the forcing or initial/boundary conditions. 

A deeper understanding is achieved by considering the Kolmogorov hypothesis |2H [2S] 
on the universality of statistics of velocity multipliers studied in the context of stationary 
turbulence. For the Sabra model, this hypothesis reduces to considering the absolute ratios 
Wn and the phases A„ dehned as [261127]. 

Wn = \Un/Un-l\, A„ = arg(M„_2Mn-lO- (9) 

The two variables z„ = {wn, A„) determine multiplicatively the shell speed for the known 
speeds Un-i and Un -2 as 

Un = Un-lWne~"^"^'-^"-\ dn-2 = argM„_2. (10) 

It is conjectured m that the probability distribution of the multiplier variables 

/^(. . . , Zn—\i Z^, Zn^\i ■ . .) (H) 

is universal in the inertial interval, which is conhrmed by numerical simulations. Further¬ 
more, the distribution 0 is short-range (correlations between z„ and Zn+j decay rapidly 
with increasing Ijj) and it is homogeneous in n (invariant under the shift z„ i—)■ Zn+j for any 
j). Our numerical simulations lead to the same conclusions, see Fig. 

The universality demonstrated in Fig. (also confirmed for the phases A„) is, in fact, 
a stronger result, because it is verihed at every time for the instantaneous non-stationary 
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FIG. 5. (a) Probability density function for the shell speed multipliers Wn = \un/un-i\- (b) 

Correlation coefficients corr(u)„, rcn+j). The results are shown at times t = 1, 1.5, 2 (green, red 
and blue) with n = 12 ,..., 15. Collapse of all the curves confirms universality of the probability 
distribution for velocity multipliers at small scales. 

probability distribution, i.e., not in the stochastic equilibrium. This means that the sponta¬ 
neously stochastic distribution of the Sabra model, written for the multipliers, has a universal 
form at the side of small scales. This can be stated in a different way: the non-stationary 
probability distribution (defined in the limit Re —)■ cxd) evolves in time under the small-scale 
boundary condition given by the universal multipliers distribution ( [TT| ). We now show how 
this boundary condition can be implemented in order reduce the problem to a purely inviscid 
evolution of the probability distribution. 


VII. STOCHASTIC SMALL-SCALE BOUNDARY CONDITION FOR INVISCID 
TURBULENT FLOW 


The evolution of a statistical distribution P„(mi, ...,-u„) for the purely inviscid Sabra 
model is governed by the continuity equation 


dPn 

dt 




k=l 


d d 

((ft^l'Ui, . . . , Un)Pn) + • • • ) '^n)Pn) 


du 


= 0 , 


( 12 ) 
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where we denoted = Rew^ and u]^ = Imn^, and ... ,Mn) means the conditional 

average of iff. for hxed Ui,... ,Un and time t. The derivatives iik = iff. + iii\. are given by 
Eq. ([^ with the vanishing viscosity z/ = 0. Since the right-hand side of Eq. ([^ depends 
on the two neighboring shell speeds on each side, one needs valnes of the speeds m_i, mq, 


Un+i and Un +2 in order to close Eq. (12). Since the speeds m_i and uq mimic the large-scale 


physical boundary, they should be deterministic (recall that, in our simulations, we chose 
M_i = Uo = 0). On the contrary, the small-scale speeds Un+i and Mn +2 are determined by 
the turbulent dynamics in the inertial interval, so they are stochastic. 

Recall that every shell speed can be given as a function Un = Un(zi,..., z„) for given uq 
and M_i with the recurrence relations (10). Thus, we can to do the closure of Eq. (IT^ in 


terms of the probability distribution ■ ■ ■, ^n) for the multiplier variables z„ = (wn, A„), 

using the exact relation 

Pn+2(^ly ■ ■ ■ 1 '^ni ^n+l) ^ 91 + 2 ) Pnif'li • • • 1 '^ri^Pcif‘n+l\'^ni 1; • • ■)Pc(.'^n+2\'^n+lj • • •)> 

(13) 

where Pc{zn\zn-i, z„_ 2 ,...) is a conditional probability. If n belongs to the inertial interval, 
the function Pc is universal, independent of n and has a rapidly vanishing dependence of 
Zn-j for increasing j, as we argued in the previous section. This means that Pc can be 
considered to be given, e.g., from numerical simulations. This provides an explicit closure. 


written as a continuity equation for similarly to Eq. (|12|) as 

dp 

- + 

k=l 


dt 


n p 

+ X] / C^^^^n+lC^^^^r^+2dA„+l(iA„+2 

7_1 ^ 


d 

dWr, 


fn.Pj 


n+2 


+ 


d 


dAr. 


A„P 


n+2 


= 0 , (14) 


where Pn +2 is given by Eq. (|13|), and all derivatives Wn and A„ are defined explicitly in 

j; note a notation difference due to the 


terms of Zi,..., z „+2 (the formulas are given in 
factor kl/^). For the probability distribution obtained in the limit Re —)■ oo, the described 
purely inviscid formulation is expected to be exact in the limit of large n. 

We arrived at a closed formulation for the evolution of a probability distribution, which 
is valid for a purely inviscid flow, in the limit Re —)• cx). Instead of the viscous regularization 


corresponding to hnite Reynolds numbers, the inviscid continuity equation (14) features 


the stochastic anomaly governed by the probabilistic “boundary condition” (13) at small 


scales (large n). Since this condition is written in terms of random multipliers, it is natural 
that the stochastic component enters the flow spontaneously after the blowup when all the 
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small-scale speeds Un get excited, see Fig. |^a). Now we are ready to verify the proposed 
stochastic dehnition of inviscid dynamics numerically. 


VIII. PURELY INVISCID TURBULENCE: NUMERICAL TESTS 


We simulate the purely inviscid evolution of the probability distribution by comput¬ 
ing numerically 10^ random sample solutions of the inviscid Sabra model for the shells 
Ui{t),..., Un{t). We consider the initial time t = 0.8 and take the initial, boundary and 
forcing conditions from the viscous simulations for Re = 10^°, which were described in Sec¬ 


tions IV and VI and will be used for comparison. Our initial time t = 0.8 is relatively close to 


the blowup and one can see in Fig. |^a) that the inertial interval is well established starting 
from the shell n = 7. Thus, we perform two different numerical tests, one with the minimum 
number of simulated shells n = 7 and the other with n = 10. In each simulation, the inte¬ 
gration is carried out using the explicit second-order Adams-Bashforth method, where the 
shells Un+iit) and Un+ 2 it) are chosen randomly at each time step according to the universal 
conditional probability Pc as we explain below. 

The universal conditional probability is implemented in terms of the multipliers = 
(wfc, Afc), which uniquely dehne the shell speeds, as explained in Section [Vn Since corre¬ 
lations between z^, decay rapidly with a shell separation, see Fig. [^b), we limit our con¬ 
sideration to the four subsequent shell numbers in Pc(zfc|zfc_i, Zfc_ 2 , We determine Pc 

numerically by using 10 ^ random samples (zi, Z 2 , Z 3 , Z 4 ) = (zfc_ 3 ,..., z^), which are taken at 
different times in the inertial interval for a single numerical simulation of the Sabra model. 
This data is computed once and stored. Then, at each time step of the integration method, 
given the multipliers z„_ 2 (t), z„_i(t), Zn(t), one choses the nearest triple zi, Z 2 , Z 3 and assigns 
the corresponding value Zn+i{t) = Z 4 . Next, the same procedure is applied to dehne z„+ 2 (t). 
These multipliers dehne the shell speeds Un+i(t) and Un+ 2 (t) necessary for computing the 
right-hand sides of the Sabra model equations (|^. 

The last column in Fig. [^presents the simulation results performed with a total number 
of n = 7 shells. The probability density of the second shell speed U 2 {t) is shown on complex 
plane at three diherent times. It is quite remarkable that with only 7 shells, i.e., ehectively 
with the shells of the forcing range only, the inviscid probabilistic formulation reproduces 
accurately the results previously obtained with very high Reynolds numbers, see the second 
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FIG. 6. (a) Mean values and (b) standard deviations for real parts of the large-scale shell speeds 
ui and U 2 as functions of time. Bold black curves correspond to the solution with Re = 10^*^ (same 
as in Fig.[^. Solid red (dashed green) curves correspond to the purely inviscid stochastic evolution 
simulated with n = 10 (n = 7) shells. The blue dotted curve in plot (b) corresponds to a different 
random choice of Un+i and Un +2 (for n = 10) highlighting the strong dependence on a stochastic 
small-scale boundary condition. 


and third columns in Fig. [Tj For quantitative comparison we plot in Fig. [^the mean values 
and standard deviations of the real parts, ReMi(t) and ReM 2 (t), as functions of time. Here 


the bold black curves correspond to numerical simulations of Section IV with the Reynolds 
number Re = 10^*^ (see also Fig. |^, which are compared with the results of the inviscid 
probabilistic formulation with n = 10 shells (solid red lines) and n = 7 shells (dashed green 
lines). One can see that the results for the mean values are almost undistinguishable. The 
standard deviations in the case n = 7 have a small difference, which decreases substantially 
for n = 10. 

We conclude that the obtained numerical results strongly support the theoretical con¬ 


struction of Section VII In order to see that the universal stochastic boundary condition rep¬ 
resents a hne “tuning” of the solution in the inertial interval, we performed simulations with 
a different boundary condition. Namely, we considered the shells Un+i{t) = 2“^/^e*®W„(f) 
and Un+ 2 {t) = which satisfy the Kolmogorov scaling and have independent 

uniformly distributed random phases 9i and 02. The resulting standard deviation of the real 
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part Reu 2 {t) is presented in Fig. [^b) by the dotted blue curve, which rapidly diverges from 
the high-Reynolds-number solution. 


IX. CONCLUSIONS 

The problem of turbulence for incompressible flows at large Reynolds numbers can be 
formulated as predicting large-scale flow through an adequate description of small-scale 
fluctuations. In this work, we showed that even if forcing, initial and boundary conditions 
are deterministic, the physically relevant flow description leads to intrinsically stochastic 
dynamics. The corresponding probability distribution is well-dehned at each finite time in 
the limit of large Reynolds numbers if one uses a regularization term with a random viscosity. 
In this way, the physically relevant stochastic solutions follow in the limit of vanishing 
(random) viscosity, similarly to discontinuous solutions for the Burgers equation that appear 
in the limit of vanishing (deterministic) viscosity. We also show that the developed turbulent 
flow allows a purely inviscid description, which is dehned for the probability distribution 
by means of a universal stochastic “boundary condition” at small scales. An important 
implication is that such a probability distribution can be predicted (computed) at each 
hnite time. 

We provided an extensive numerical evidence that the above description of the turbu¬ 
lent flow is valid for the Sabra shell model of turbulence. We argue that the stochastic 
anomaly is fundamentally different from the hnite-dimensional chaos, where the evolution is 
intrinsically deterministic at any hnite time. On the other hand, it is closely related to the 
phenomenon of spontaneous stochasticity due to non-uniqueness of Lagrangian trajectories 
in rough (singular) velocity helds. As we demonstrated for the Sabra model, the singular 
velocity appears naturally in the evolution as a result of a hnite-time blowup, immediately 
triggering the stochastic behavior. This provides an example of a classical system whose 
hnite-time behavior can only be described with the probabilistic terminology, a property, 
which is normally attributed exclusively to quantum systems. 

Our theoretical considerations are based the Kolmogorov-Obukhov theory and the con¬ 
cept of inverse cascade of stochasticity discussed in the literature for long time mm- Since 
these ideas are valid both for the Sabra model and for the 3D Navier-Stokes equations, we 
conjecture that a similar formalism applies to the realistic turbulence for high Reynolds num- 
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bers. Though the existence of a hnite-time blowup (that may “turn on” the stochasticity) 
remains an open problem for the inviscid flow (3D Euler equations) |2H], numerical simu¬ 
lations suggest the blowup at a physical boundary [22], while nearly exponential vorticity 
growth is typical far from the boundary, e.g., [301 EH- The spontaneously stochastic distri¬ 
bution for the Navier-Stokes equations can be accessed numerically by performing a large 
number of simulations for the same deterministic large-scale conditions and small random 
viscosities, which is a challenging problem due to very high requirements for computational 
resources. 
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